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Abstract 

We study non-Gaussian properties of the isocurvature perturbations in the dark 
radiation, which consists of the active neutrinos and extra hght species, if exist. 
We first derive expressions for the bispectra of primordial perturbations which are 
mixtures of curvature and dark radiation isocurvature perturbations. We also discuss 
CMB bispectra produced in our model and forecast CMB constraints on the non- 
linearity parameters based on the Fisher matrix analysis. Some concrete particle 
physics motivated models are presented in which large isocurvature perturbations 
in extra light species and/or the neutrino density isocurvature perturbations as well 
as their non-Gaussianities may be generated. Thus detections of non-Gaussianity in 
the dark radiation isocurvature perturbation will give us an opportunity to identify 
the origin of extra light species and lepton asymmetry. 



1 Introduction 



Recently, several cosmological observations independently suggest that the effective num- 
ber of neutrino species in the Universe is larger than the standard value, i.e. ANes = 
— 3.04 ^ 1. According to recent observations of the primordial abundances of light 
elements, it is constrained as A'cfr = S.GStgyQ at 2a level (with slight dependence on the 
center value on the measured neutron lifetime) [1]. On the other hand, recent observa- 
tions of the cosmic microwave background (CMB) anisotropy at small scales in combina- 
tion with WMAP [2] and standard distance rulers [3, 4, 5], give A^efr = 4.56 ± 0.75 [6] 
and A^eff = 3.86 ± 0.42 [7] at Icr level. These results may be evidences for the exis- 
tence of extra radiation component, other than the three species of active neutrinos, 
in the Universe. See also Refs. [8, 9] for limits on the mass of extra radiation compo- 
nent. Motivated by these observations, models for explaining AAefr — 1 were proposed 
[10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. The Planck and other projected CMB obser- 
vations will improve constraints on N^s by an order of magnitude (see e.g. Refs. [22, 23]), 
and AAcff — 1 can be clearly tested in the near future. 

Once it will be proven that extra radiation indeed exists, it is important to under- 
stand the origin of extra radiation in the early Universe. In the previous work [17], it is 
argued that one way to probe this is to see how they fluctuate at the cosmological scales. 
Observationally, the extra radiation and neutrinos are not discriminable and we call the 
mixed fluid of them as dark radiation (DR) . The DR can have isocurvature perturbations, 
depending on how they are produced in the inflationary Universe. It was shown that the 
DR isocurvature perturbations affect the CMB anisotropy and constraints on the ampli- 
tude of isocurvature perturbations were derived using recent CMB and other cosmological 
observations. Future forecasts on the constraint were discussed in Ref. [24]. 

While perturbations are assumed to be Gaussian in the most part of Ref. [17], the 
possibility of large non-Gaussianities in the DR isocurvature perturbations was also briefly 
pointed out. In this paper, we present detailed analysis on these non-Gaussianities. Non- 
Gaussianities in the DR isocurvature modes would have rich information on the properties 
of DR. We note that there are several studies on non-Gaussianities in the cold dark matter 
(CDM) and baryon isocurvature perturbations [25, 26, 27, 28, 29, 30, 31, 32, 33, 34]. 
However, this is the first paper that studies non-Gaussianities in the DR isocurvature 
perturbations, including those in the neutrino density isocurvature perturbations. In this 
paper we focus on the local type non-Gaussianities at bispectrum level. 

The paper is organized as follows: We first present bispectrum generated from mixtures 
of primordial DR isocurvature and adiabatic perturbations in Section 2. In Section 3, we 
apply these results to the CMB angular bispectrum and discuss how non-Gaussianities 
in DR isocurvature perturbations manifest in the CMB anisotropy. Then we discuss 
constraints on these non-Gaussianities from CMB observations in Section 4. Here we 
forecast constraints from the Planck satellite and an ideal survey limited by the cosmic 
variance, based on the Fisher matrix analysis. We mention some particle physics models 
which may lead to large isocurvature perturbations in the extra radiation and neutrinos as 
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well as non-Gaussianities in them in Section 5. The final section is devoted to summary. 



2 Non-Gaussian curvature and isocurvature pertur- 
bations 

In this section we derive formulae for the (non-Gaussian) curvature/isocurvature pertur- 
bations based on the 5iV-formalism [35, 36]. We mostly follow formalsm in Ref. [17]. 
Various modes of primordial perturbations, including the adiabatic mode ( and some kind 
of isocurvature mode S, can be generated from fluctuations in scalar flelds and given as 

C = N^Mi + ^N^^^^S(j)Mj + (1) 

S = S^^6(f)i + ^S^^^^6(f)i6(f)j + — 

Here, 5(j)i is quantum fluctuation of a scalar fleld (pi, whose mass is smaller than the 
Hubble parameter during inflation, H[^f. Hereafter, we concentrate on the isocurvature 
perturbation in the dark radiation (DR) denoted by S'dr. We here again emphasize that 
the DR consists of both active neutrinos and extra light particle species. 

We can express the power spectra of the auto- and cross- correlation functions of ( and 
S'dr as follows, 

(C(^i)C(4)) = (27r)=^5(^i + 4)P«(A;i), 
(C(^i)5dr(4)) = i2nf5{h + k2)P'^'^Hki), (2) 
(^dr(^i)5dr(4)) = {2nf5{h + 4)P^°«^°«(A;i). 

where 

P^'^m = N^^S^^Ps^ik), (3) 
pWoa(^) = S^Ps^k), 

Here, we have neglected higher order terms and P5^{k) is the power spectrum of the 
fluctuations of the scalar flelds, 

{5Uki)H3{k2)) = {27rf6{h + k2)P5^{ki)6,j, (4) 

where is the scalar spectral index*^ and ko is the pivot scale chosen as fco = 0.002Mpc"^. 
Note that the above power spectra and correlation have same spectral shape up to this 



The scalar spectral indices for (f> and a do not coincide in general. In the following we assume they 
are the same just for simplicity. 
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order. We therefore adopt = P'"'' as the normahzation of power spectra and other 
power spectra can be expressed in the form P^ times some constants hereafter. 
The bispectra of ( and S'dr are defined by the following equations, 



(C(fci)C(fc2)C(4)) = (27r)35(fci + k2 + h)B^aih, k2, h), 
(C(^i)C(^2)^dr(4)) = {2Txf5{h + k2 + h)Bccs{h, h, h), 
(C(^i)5dr(4)5dr(4)) = {27r)'Sih + k2 + h)Bcss{ki, k^, ks), 
{SBRih)Sr,Kih)SMh)) = {2Ti)H{k^ + k2 + h)Bsss{ku k^, k,), 



(6) 



so that the primordial bispectrum can be written in the form of 

P^^^^^«(fci, k2, k,) = f^f'^'^'ik,, k2, ks)Pdk2)Pdh) + (2 cyclics of {123}), (7) 

where the each of the subscript Ai {i = 1, 2, 3) is either ( or S'dr- The coefficients 
/nl'"^^^^ represent magnitudes of non-Gaussianities and in the following we call them non- 
Gaussianity parameters. Note that our definition of f^l'^^^'^ is consistent with Ref. [34] 
besides difference in types of isocurvature perturbations considered. We also note that 
if there is only a single scalar field which sources primordial perturbations and there are 
only adiabatic perturbations, /^l^ is related to the ordinary non-Gaussianity parameter 
/nl via 

/^f = IhL. (8) 
By using the expansion (1), we can explicitly write the non-Gaussianity parameters as 

/nl - ^Nlf + ^CM^^^)' 

/nl - /nl - tf,j2 \2 ^ (AT2 mi^/CfelvJ , 



/nl = — T^rzYi — + imV^ ^^y.ki,L), 



(9) 



/nl - /nl - (^2^)2 + ^^ik,L), 

Q Q Q Q Q Q 

where = {k^ /27r^)Pt^{k) is the dimensionless power spectrum of the curvature pertur- 
bation, kb = min{/ci, k2, k^} and L is the infrared cutoff scale [37, 38], which should be set 
to be a scale comparable to the present horizon scale. 
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3 CMB bispectrum 

CMB bispectrum from non-Gaussian curvature and extra radiation-isocurvature pertur- 
bations are to be discussed. For CDM isocurvature perturbation, similar analysis is done 
in Refs. [25, 28, 29, 31, 34], which extend the analysis of Ref. [40] to isocurvature per- 
turbations. While we consider the isocurvature perturbations in extra or dark radiation 
'S'dr; we include not only temperature anisotropy but also polarization one, which is not 
included in these works (see also e.g. Ref. [41]). 

First, we denote primordial perturbations by X^, where the subscript A is either ^ or 
5'dr. CMB anisotropy is given by 



''Im 



r ^3 J, 

4^(-^)' E J j^39rimi{k)x^, (10) 



where the subscript P represents the type of CMB anisotropy and should be either T or 
E, and gf^i^k) is the transfer function at linear order. 

First the power spectrum of the CMB anisotropy is expressed as 

= {a^>f^>). (11) 
It is given in terms of the primordial perturbations as 

C^^P^ = -J2[ k''dkgf^'''{k)gf'''\k)P^^^^{k). (12) 
^ A1A2 

where P'^^^'^{k) is the power spectrum of in the wave number space defined in Eq. (2), 
which is conveniently written as 

(X^n^i)^*^'(^2)) = P^^^HA;i)(27r)35(3)(^^ _ k^). (13) 
Let us now consider the bisepctrum of CMB anisotropy in the harmonic space, 

^PiPaPs ^ / Pi P2 P3 \ nA\ 

Using Eq. (10), B^^^f^^^^^^^ can be written as 

i=i L J y^"") 

x5^i^^^«(A;i, k2, k,){27rf6^^\h + h + h), (15) 

where B^^^'^^^{ki,k2,k3) is the bispectrum of X^ in the wave number space defined in 
Eq. (7), which are conveniently written as 

{X^^{h)X^^{h)X^^{h)) = B^^^^^-^{kuk2,ks){2nf6'''\h + ^2 + h)- (16) 
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Due to the statistical isotropy assumed here, _B^i^2^3 jg independent of ki 
Using the following formulae 



dkY;jk)e^'-''' = Am^ji{kr)Y;jf), 



(17) 
(18) 



-'^Zimi/2m2'3'Ti3 



- / k'tdhg^^''^{h)jiXhr) 



Eqs. (15) can be reduced into 

r ^ 

xB^'^^^^^ik^k^ks) I drUi'^CmS^^] 
In Eq. (19), we can factor out the Gaunt integral 



1=1 



(19) 



i=l 



(20) 



which manifests the statistical isotropy. In terms of the Wigner-3j symbol, Qi^l^i^ ^ can 
be rewritten as 



hhh 

Then we obtain 



(2/i + l)(2/2 + l)(2/3+l) fh k l3\ ( h k h 

A-K I / I mi m2 msj ' 



nPiP2Ps _ ^m,im2m3 J P1P2P3 

-'-'limi, 121712, hm^ ^khh hhh ' 



where b^^i^f^^ is the reduced bispectrum given by 

3 



A1A2A3 ■ 



B 



A1A2A3 



(21) 
(22) 

ik,,k2,k3) (23) 



This is the most general expression of CMB bispectrum, in the presence of non-adiabatic 
primordial scalar perturbations. This is applicable to any types of primordial non-Gaussianities. 

Now we focus on primordial perturbations with local-type non-Gaussianity, which can 
be written in the form of Eq. (1). 

Given the primordial bispectrum of Eq. (7), the reduced CMB bispectrum b^^i^f^^ can 
be written as 



hPlP2P3 

%hh 



A2A3^A.P,A2P2A3P3 ^ ^2 cyclics of {123}) 



(24) 



A1A2A3 



Here, b^]^,^^'^'^'^'^' is defined as 



where 



^AlPl,A2P2A3P3 



AP 



a, r 



^ I edkgfP{k)ji{kr) 



APi 



(25) 

(26) 
(27) 



As long as only two kinds of initial perturbations ( and Sdr are included, there are 
only six independent non-Gaussian parameters [34]. For latter convenience, we denote 



them by f^l {a 



6), which is defined as 



.(1) 



/nL — JNL ' 



f(2) 



■Sdr,CC 



/nl 



f(3) 



(5) 



f ■S'driC'S'dr /S'driS'drC f 



(6) 



We also pile up bf^i^l^^'^^'^'^^^^ into six types of reduced bispectra ftj^j^zg 



/C.C'S'dR /Cj'^drC 

iNL ~ iNL ' 

f S'drjS'drS'dr 

(a) P1P2P3. 



^(1) P1P2P3 

^(2) P1P2P3 
'l'2'3 

^(3) P1P2P3 

'l'2'3 
^(4) P1P2P3 

,(5) P1P2P3 
"hhh 

^(6) P1P2P3 

'l'2'3 



C2r'^'^' + (2 cyclics of {123}), 

+ (2 cychcs of {123}), 
+ + (2 cychcs of {123}), 

Ca^^''^'^^''^' + (2 cychcs of {123}), 



5drPi,C^'25drP3 

«2«3 



,SDRPl,5DRP2Cf"3 
'^hi2«3 



(2 cyclics of {123}), 



^3^'"^°''^'^°''^' + (2 cychcs of {123}) 
Then the total CMB bispectrum in Eq. (24) can be rewritten as 



(28) 
(29) 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



LP1P2P3 _ \^ fWtW 1 



(a), (a) P1P2P3 



a=l 



Fig. 1 shows the temperature bispectra b['^}'^i^^ in isosceles triangular configurations 
with li = h- Cosmo logical parameters adopted here are the mean parameters for the fiat 
power-law ACDM model from the WMAP 7-year result [2]. In numerical calculation, the 
transfer functions gf^{k) are computed using the CAMB code [42]. Among six bispectra 



(36) 



(a) 



and b 



(3) 



'hhh^ "hhh "^^^"^ "^1^2/3 ^ends to be larger than others in most configurations, although 
configurations shown in the figure are limited. On the other hand, b[^j^i,^ is in general 
the smallest. We also see in the figure that the acoustic peaks and troughs in different 




Figure 1: Shown are the temperature bispectra b^ij^i^'^ ■ ^txhi^^'^ (sohd red), hfj^j^^'^ (short- 
dashed green), hf^^^'^ (dotted blue), hf^^i^'^ (dot-dashed), ^if^^a^^ (long-dashed), ^Ss^^ 
(dot-dot-dashed) in isosceles triangular configurations with li = I2 are plotted as function 
of li with fixed l^. I3 is set to 10 (top left), 50 (top middle), 100 (top right), 500 (bottom 
left), 1000 (bottom middle), 1500 (bottom right). In each panel, the shaded region at low 
multipoles shows configurations which fail to satisfy the triangular condition i.e. |/i — /2I > 
h>h + k- 
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Figure 2: Same figure as in Fig. 1 but for ^["/j^J^- 
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Figure 3: Same figure as in Fig. 1 but for ^^ij^if^ ■ 
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bispectra show up at similar multipoles. This is because both in the adiabatic and neutrino 
density isocurvature (NID) modes, the acoustic oscillation of the photon-baryon fluid has 
cosine-like phase (For detailed discussion on the acoustic oscillation in the NID mode, we 
refer to Ref. [17]). This is contrastive compared with the bispectra from non-Gaussianities 
in the matter isocurvature modes (See Ref. [28]). While oscillating features in the different 
bispectra are more or less similar, the global spectral shapes are different, which allows us 
to distinguish one from another. 

We also plotted the bispectra arising from the correlation of two temperature and 
one polarization anisotropies, b^'^j^ij^ and b^fi'^i^'^ in Figs. 2 and 3, respectively. From 
these figures, we can see that above discussions on spectral shape of the bispectra are still 
qualitatively true when polarization is included. 



4 Forecast for a CMB constraint 

As we have seen, CMB bispectrum arising from the non-Gaussian isocurvature perturba- 
tions in extra radiation is distinct from the usual one from the non-Gaussian curvature 
perturbations. Therefore we can discriminate different non-Gaussinities in primordial per- 
turbations from the observation of CMB anisotropy. To discuss this issue in a quantitative 
manner, we perform a Fisher matrix analysis. 

In the limit of weak non-Gaussinity, the Fisher matrix for the non-Gaussianity param- 
eters fl^l is given by [40, 41, 43] 



h<l2<h 

X E E 

P1P2P3 Q1Q2Q3 



(2/i + l)(2/2 + l)(2/3 + l) 



An 



h h h 




,(a) PlP2P3tn^^~l-\PlP2P3\QlQ2Q3,{b) Q1Q2Q3 

\i2i3 y^^^ nihh \i2i3 



(37) 
(38) 



where [Gov" 



is the inverse covariance matrix. Assuming that the ob- 



served sky coverage is unity and the instrumental noise is isotropic, the covariance matrix 
[Cov]^^]^^'^^'^''^' can be given as 



[cov]sr' 



1Q2 



A 



pPlQlpP2Q2p 



P3Q3 
I3 ' 



(39) 



where C/^*^ = Cf'^ + A*";^*^ is the total angular power spectrum, which is the sum of ones 
from the CMB C/''^ and instrumental noise A^^*^. A/j/j^g takes values 6, 2, 1 for the cases 
that all /'s are the same, only two of them are the same and otherwise, respectively. 
Following [44], the noise power spectrum A^^^*^ can be approximated as 



A^, 



PQ 



^PQ^FWHmC^P 6Xp 



/(/ + 1 



fl2 

I^FWHM 

81n2 



(40) 
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Table 1: Expected uncertainties for non-Gaussianity parameters A/j^^ the case with 



survey 


Jnl 


f(2) 


f(3) 


f(4) 


f(5) 


f(6) 


Planck 


22 


101 


21 


116 


163 


164 


CVL 


3.5 


14.0 


3.7 


15.9 


15.4 


17.3 



where ^fwhm is the full width at half maximum of the Gaussian beam, and crp is the root 
mean square of the instrumental noise par pixel. For cases of multi-frequency observations, 
is given via the quadrature sum over all the frequency bands. 
In the analysis, we consider only CMB bispectra from primordial non-Gaussianities, 
assuming contaminations from other sources are negligible. Since many of these sources 
including point sources and the thermal Sunyaev-Zel'dovich (SZ) effect have frequency 
spectra different from the black body, above assumption can be to some extent achieved 
by exploiting observations at multi-frequency bands. Other contaminations such as the 
lensing of CMB, the kinetic SZ effect, and the patchy reionization would not affect our 
results significantly. 

As our models have six non-Gaussianity parameters fl^^, the full Fisher matrix Fab is 
a 6 X 6 matrix. However, it may sometimes occur that some of the parameters are not of 
primary interest and we want them to be marginalized over. In such the case, the Fisher 
matrix for the remaining non-Gaussianity parameters can be given as the inverse of the 
principal sub- matrix of the inverted full Fisher matrix [45]. 

In Figs. 4 and 5, shown are 2-dimensional constraints on the non-Gaussianity parame- 
ters expected for Planck and a cosmic variance limited (CVL) surveys. Here we fixed N^, 
to 4, which is suggested by recent observations we mentioned in Introduction. On each 
panel, constraints on a pair of f^^^ are shown; other four non-Gaussianity parameters are 
marginalized over in Fig. 4 while they are fixed to zero in Fig. 5. Hereafter we will refer 
to constraints shown in Fig. 4 and 5 as marginalized and non-marginalized constraints, 
respectively. 

In Table 1, we listed expected uncertainties in the non-Gaussian parameters A/jlj^, 
which are defined by 

A/S = F-a'- (41) 
From the table as well as figures, we can see that among the six non-Gaussian parameters 
/nl' /nl ^^"i /nl h^ constrained tighter than others. Planck (a CVL survey) can 
constrain f^l and f^l to about 20 (4). On the other hand, expected constraints on other 
are weaker with factor from five or eight. This result is consistent with our discussion 
in the previous section, where we showed that in the squeezed configurations, i>i^i'^if^^^ 

and bf^i^i^^^^^ are in general larger than other bispecra. We can also see that a CVL survey 
can significantly improve the constrains on all non-Gaussianity parameters from Planck 
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Table 2: Same as in Fig. 1, but for the case with Ncs = 3.04. 



survey 


Jnl 


f(2) 
JNL 


f(3) 
JNL 


f(4) 
JNL 


.(5) 
JNL 


/■(6) 


Planck 


21 


126 


27 


187 


257 


339 


CVL 


3.5 


18.3 


5.0 


27.2 


26.4 


39.3 



by an order of magnitude. On the other hand, in the case of the matter isocurvature 
mode, constraints on some of non-Gaussianity parameters improve little as we measure 
higher and higher multipoles as shown in Ref. [34]. This difference reflects that CMB 
anisotropics at high multipoles are damped in the case of the matter isocurvature mode, 
while they are comparable in amplitude with the adiabatic mode in the extra radiation 
isocurvature mode. 

We also performed the same analysis for the case of fixed = 3.04 and marginalized 
and non-marginalized constraints on are shown in Figs. 6 and 7, respectively. Param- 
eter uncertainties A/j^j^ are listed in Table 2. In the context of isocurvature perturbations 
in extra radiation, such the case is realized when, while the fraction of extra radiation 
in energy density of DR is quite small, the amplitude of isocurvature perturbations in 
extra radiation Sx is large enough for the total DR isocurvature perturbations Sbr to 
be yet non- negligible. On the other hand, this is also naturally realized without extra 
radiation; S'dr is non-zero if there are isocurvature perturbations in the lepton number 
and non-Gaussian isocurvature perturbations in the lepton number may be produced in 
the Affleck-Dine mechanism, as shown in Sec. 5.2. Compared with the case of A^j, = 4, 
constraints on f^l are less stringent for the case of N^, = 3.04. This is simply because ef- 
fects of DR isocurvature perturbations becomes less significant as A*";, becomes small given 
fixed S'dr. However, there are qualitatively little difference between the cases with = 4 
and 3.04 and parameter uncertainties increases by factor 2 at most. 

5 Models for non-Gaussian isocurvature perturbations 
in dark radiation 

In this section, we refer to some of particle physics models in which the DR isocurvature 
perturbations and their non-Gaussianities arise. We discuss two cases separately. In one 
scenario, the DR isocurvature perturbation is carried by extra light species. In the other 
scenario, ordinary neutrinos have large isocurvature perturbation. 

5.1 Extra light species 

Let us consider the cosmological scenario considered in Ref. [17] where two scalars, the 
inflaton and the curvaton cr, which is light during inflation, contribute to both the 
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Figure 4: 2d marginalized constraints on the non-Gaussianity parameters /j^^ expected 
for Planck (solid red) and CVL (dashed green) surveys. is fixed to 4. Inner and outer 
contours correspond to constraints at 1 and 2 a levels. 
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Figure 5: Same figure as in Fig. 4, but the non-marginalized constraints are shown here. 
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Figure 7: 2d non-marginalized constraints on non-Gaussianity parameters for = 3.04. 
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adiabatic and isocurvature perturbations. Using the SN formalism [35, 36], the expansion 
coefficients of (, such as A'',^, A^'^., are given as 



1 V 
3 + R f RrR^^^ 



6(Ti 

1 

Ml 



1 - 



vv,, 



+ 



rixJ^x 



R 



X 



Vi 



2 3 + i? / RrRr ^ , RxRx^ 



9at 4 



Rr 



+ 



R 



X 



3 + 4i?-2i?2 -2(3 + i?) 



In a similar manner, those of S'dr are given as 



Rr-R 



Rr 



(42) 



KxrCx 



R 



X 



3 -\- R RrRx I Rr 



2ai R 



DR 



DR 



Rr 



3 -\- R RrRx I R. 



(^) 



Rr 



X 




R 



DR 



Rr 



0. 



(43) 



RxRx\^ + -Rdr) 
Rx 



The meanings of the symbols are as follows: Mp is the reduced Planck mass. V is the 
potential of (p. and V^^ are the first and second derivatives of V , respectively. Ri 
is the ratio of the energy density of a fluid i to the total energy density at the decay 
of cr, and Ri is that at the electron-positron annihilation. The subscripts r, X and DR 
mean the relativistic particles in the Standard Model, the extra radiation and the dark 
radiation, respectively. R^^^ is the ratio of energy density of the fluid i generated by a 
decay at that time, and R = 3Ra/{4: — Ra), where R^ is the ratio of energy density of a 
to the total energy density at its decay, cij is the amplitude of the oscillation of cr when 
it starts to oscillate, ^ 0.405 is the ratio of the energy density of neutrino to that of 
standard model relativistic particles (photons and neutrinos) after the electron-positron 
annihilation. Using these quantities, we obtain 



AN, 



3R 



X 



cS 



(44) 
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We refer to Ref. [17] for details and derivations of these quantities. 

Using these quantities, the non-Gaussianity parameters defined in Eq. (9) are expressed 

as 

f^iHh, h, k,) ^ ^*ffi;^°r + Hk.L)Al(k,), (45) 



C ]\t2 q 



h, h) = f'^'^^ih, h, h) (46) 



fS^'^'^^'iku h) = f-^'^-^^ih, ^2, h) (47) 
- ln(/cbL)A.(/i;i), 



C2 C C3 

(ivJTW MTiv,^) 



Here, we include the "quadratic" type components [39, 25, 28, 29, 31], which consist 
of three quadratic terms of 50j in each ( or S'dr, in addition to the leading "linear" 
components, which consist of a quadratic term in one of three ( or 5'dr and two linear terms 
from others. While /^l^, etc. in Eqs. (45)-(48) are in principle not constant, their scale- 
dependences due from the factor ln(A;bL)A^(/cj) are quite moderate since P({k) is nearly 
scale-invariant. Therefore, so long as we consider observations sensitive to scales over only 
a few orders of magnitude, we can approximately ignore the scale-dependences. When we 
in the next section consider CMB signatures of non-Gaussian isocurvature perturbations 
in dark radiation, we adopt this approximation and regard the quantities f^'^'^, etc. as 
constants. Then given constant f^'['", etc., as shown in Eq. (7), the bispectra can be 
written as 

B<<<{ki, k2, ks) = /^f [Pc(fc2)Pc(^3) + (2 cyclics of {123})] , (49) 
B^^'^^^ik^, k2, ks) = /^f °^ [Pcik2)Pdh) + Pdks)Pcik,)] + /^r'^^n(^i)n(^2), (50) 

5C5Da5„H(^^^ ^2, ks) = /^f°«^°«Pc(^2)Pc(^3) (51) 

+ fSun,Sun( [p^^ks)Pdh) + ^(^1)^(^2)] , 
55dr5dr5dr(^^^ f.^^ _ j5dr,5dr5dr [p^(^2)Pc(^3) + (2 CyclicS of {123})] . (52) 

In the following subsections we consider two cases. One is the case where the a domi- 
nantly decays into extra light species X, while the curvature perturbation is dominantly 
generated by the inflaton (Sec. 5.1.1). The other is the case where the a decays into 
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ordinary radiation and is the dominant source of the adiabatic perturbation, while extra 
hght species are produced in thermal bath after the inflaton decay (Sec. 5.1.2). Both cases 
are realized in the framework of supersymmetric (SUSY) axion model [46] as mentioned 
in the previous work [17], and originally in Ref. [10]. We shall partly repeat discussions 
there. 



5.1.1 Dark radiation from particle decay 

Let us assume that the primordial curvature perturbation is dominantly produced by the 
inflaton : A^^ ^ iVo-,*^ and the inflaton decay only to the visible sector. This makes the 
isocurvature mode uncorrelated with the adiabatic mode. In this setup, we can approxi- 
mate parameters as R^"^ ~ 0, and Rx = R^x — Ra — 4i?/3. We also assume R^ < I since 
otherwise a dominates the Universe before it decays and the Universe would be dominated 
by X. Then Eqs. (42) and (43) are simplified as 

^ ^Rx, 

1 - 
N ~ Rv 

* (7(7 2 ^ ' 

r; _ " 1 _ ^ (53) 

2af Cy Ci, 

Here and hereafter, we assume that the inflaton does not induce the non-Gaussianity, that 
is, N^^ ^ 0. The quantity Rx is related to Rx as [17] 



1/3 

Rx. (54) 



Here 5 = ((11/4)^/^ — 1)(1 — c,^), Ci, = Pu/Pr = 21/43 is the ratio of the energy of neutrinos 
to that of all visible matters at neutrino decoupling, g*{H = T^) is the relativistic degrees 
of freedom at that time and g*{H = T^) is that at a decay. From Eq. (44), we can also 
express the effective number of neutrino species as 



'cff 



6) yg^H = r. 



1/3 



R. (55) 



As shown in (53), is comparable to Ncr in this model, then A^^ ^ is required in order to avoid 
the isocurvature mode comparable to the adiabatic mode. 
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Thus ANcs ~ 1 if -R is not much smaller than one. Using (53), we get the relations among 
the non-Gaussianity parameters as 



JnL — JNL — - JnL ' 

fCf^nS^n ^ fSun,CSun ^ g (^^^^ /^f , (56) 

/nl -^'y J /nl • 

Thus these non-Gaussianity parameters are comparable. The magnitude of them is roughly 
given by 



.= (^^)*fli+«'(^)°fliA^(i,), (57) 

where e = |M|(V0/y)^ is the slow-roll parameter. It is easily found that the first term is 
of the order of (Psj^^/P^)"^ x (l/Rx) while the second term is of the order of {Psok/^cY ^ 
(A^/i?^). Therefore, the non-linearity parameter can be large enough to be probed for 
PsuH ~ -^c Thus even if Rx is very small and there is no significant deviation 

from Ncs = 3.046, the DR isocurvature mode and its non-Gaussianity may be detected. 

Now let us estimate Rx ~ R and Sbr in the SUSY KSVZ axion model [48]. In a SUSY 
axion model [47], the saxion a, the scalar partner of the PQ axion, exists and has a mass 
777.0- which ranges from 0{keV) to (9(TeV), in accordance with the SUSY breaking scale. 
The saxion can have a large initial amplitude (Tj during inflation, and may obtain quantum 
fluctuations Sa ~ ifinf/27r if it is much lighter than the Hubble parameter during inflation. 
In this model, the dominant decay channel of the saxion is typically that into two axions. 
Relativistic axions produced by the saxion decay behave as an extra radiation X, since 
they are decoupled from ordinary matter almost completely. Here the R ~ 3Rcr/4: is given 
by [17] 



for <^ 1 and Tn^r > and 



for R and nicr < T^, where fa is the PQ symmetry breaking scale, Tr is the reheating 
temperature after the inflation and F,^ is the decay rate of the inflaton. Correspondingly, 
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the magnitude of the DR isocurvature perturbation is given by 



for i? <C 1 and > and 

for i? ^ 1 and mg. < F^. Here S'dr is regarded as S'dr = a/ (A;'^/27r2)P5pjj5^j^ which should 
be compared with ^JJWJ2tt^^)P^(^ ~ 5 x 10~^. Thus the magnitude of the DR isocurvature 
perturbation can be sizable. 



5.1.2 Dark radiation from thermal bath 

Next, we consider the case a takes a role of the curvaton, and hence it dominantly sources 
the adiabatic perturbation : N^- ^ N^f,. Moreover, we assume that the inflaton decays into 
ordinary radiation with a branching ratio r^, and into X with a branching ratio 1 — r^. The 
curvaton a is assumed to decay only into ordinary radiation. In this model, we make use 
of following approximations : R^^^ ~ 4i?/3), R^^ — 0- We also assume R^j <1 since 

otherwise the a decay releases huge amount of entropy and it dilutes the X abundance 
significantly. Under these assumptions, Eqs. (42) and (43) are simplified as 

~ 

?-c ^ 2R 1-c ^ ^62) 

1 - £, 2i? Q 1 - D Ar 

Cv Saf Cp 
The relation between Rx and Rx is given by Eq. (54) where Rx is given by 

Rx = il-r^)il-R.). (63) 

In this model, AN^s is given by 

3 fg,{H = r,)^'^'' 



(1-^^)- (64) 
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The relationships among the non-Gaussianity parameters become 

JNL — JNL ~ p. -fi-XjNL ' 

fCSunS... ^ fSun,CS^n ^ g (^—^\ ^ ^^/«^ (65) 



^5.r,5.r5.r ^ ('LAV^^/ac 



In this case, 

J^dr.SdrSdr ^ jCSdrSdr ^ jS'dr.CSdr ^ jSdr,CC ^ jC^drC ^ /''''''' (66) 

This is because a, which is the origin of the non-Gaussianity, dominantly decays to visible 
particles. Since we assume that the primordial curvature perturbation is dominantly 
produced by cr. the non-linearity parameter for the adiabatic perturbation is given by 

f^i^ih, h, A:3) ^ ^ + ^ \n{hL)Al{h). (67) 

We see that the non-Gaussianity in the adiabatic perturbation becomes large for R <^ 1 
while that of the DR isocurvature mode, /nl^ % is at most order unity unless Rx is close 
to unity. 

The above situation is actually realized in the SUSY DFSZ axion model [49] once 
the saxion is identified as the curvaton a and the axion as the extra light species X. In 
this model, the dominant decay channel of the saxion may be that into a Higgs boson 
pair. Therefore, the energy of saxion is almost converted to visible particles. On the 
other hand, there may be axions produced from the thermal bath during reheating. Let 
us suppose that the reheating temperature is high so that it satisfies Tr > Td, where 
Td — 10''GeV(/a/10^°GeV)^'^^^ is the temperature at the axion decoupling from thermal 
bath [50]. Thus axions are thermalized after the infiaton decay. The ratio of the axion 
energy density to the total energy density at the axion decoupling is given by g^{T = Td)~^. 
The infiaton decay branching ratio into X, 1 — r^, is replaced by the ratio of the energy 
density of axions to that of the whole radiation originating from the infiaton at the epoch 
of saxion decay. Thus it is estimated to be 

l-r, = ^^(»4^)'". (68) 
g^{T = Td) \g*{T = Td) / 

where Fg- is the saxion decay rate. From this, we see that 1 —r^ is typically much smaller 
than 1. Eventually, the amount of axions is smaller than that of visible particles whether 
they are produced by the infiaton or the saxion. Then we have Rx <^ 1 from Eq. (63). 
The ratio of the saxion energy density to the total energy density at the epoch of saxion 
decay is given by 

2 
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for rricr > and 



for < r^, where fi denotes the higgsino mass. The magnitude of the DR isocurvature 
perturbation is given by 



^DR =^ 6 X 10 



Vl06GeVy VlGeV/ \ 1^ J \W^^GeV J \10^^GeV J \fa J ' 

(71) 

for m„ > and 



fa \ f -f^inf \ f 



(72) 



for < r^, where we have used g*{T = Td) = 228.75. Thus the DR isocurvature 
perturbation can be sizable for some parameter choices even if Rx is much smaller then 
unity. 



5.2 Large lepton asymmetry 

Although the neutrinos are in thermal equilibrium before the decoupling at T ~ 1 MeV, 
the lepton number is conserved well after the electroweak symmetry breaking (EWSB) 
since the sphaleron effect [51] is suppressed. Therefore, if the asymmetry in the lepton 
(or neutrino) number, n^, — rip, is created after the EWSB, it survives thereafter.*^ This 
opens up a possibility that the (non-Gaussian) isocurvature perturbation in the lepton 
asymmetry is created after the EWSB [52]. If the lepton asymmetry is large enough, it 
contributes to the A^cfr as well as S'dr- A concrete example was given in Ref. [53], where 
it was shown that the late decay of Q-balls can create large lepton asymmetry. It is 
interesting because we do not need an extra radiation particle X to produce significant 
amount of DR isocurvature perturbation. Let us follow the arguments of Ref. [53] and 
estimate AA^'cfj and S'dr in this model. 

A large lepton asymmetry is created through the Affleck-Dine (AD) mechanism [54, 55]. 
Specifically, we make use of the LLe fiat direction (also called as the AD field) [55, 56]. It 
does not have the baryon number, and hence it can create large lepton asymmetry without 
producing too much baryon asymmetry. If the AD field fragments into Q-balls [57] in 
which almost all the lepton number is confined [58, 59, 60, 61] and they evaporate after 
the EWSB, a large lepton asymmetry is released and it does not washed out. Since the AD 
field may obtain quantum fiuctuations in its angular component during inflation [62, 63], 

The asymmetry in the charged lepton sector must be same as that in the baryon sector because of 
the electric charge conservation. We use the conventional terminology "lepton asymmetry" hereafter, but 
it actually means the asymmetry in the neutrino sector. 
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it results in the isocurvature fluctuation in the lepton asymmetry, i.e., the neutrino density 
(non-Gaussian) isocurvature perturbation. Non-Gaussianity in the baryonic isocurvature 
perturbation generated through the AD mechanism was studied in Ref. [26]. 

We denote by if) the AD field along the LLe fiat direction. It is lifted by the dimen- 
sion six operator in the superpotential, W = iIj^/{6M^) with M being the cutoff scale. 
We assume the gauge-mediated SUSY breaking model [64] in the following. The scalar 
potential for the AD field is given by [65] 

V={ml^,-cH'M' + M'^(^\og^^ +[^mmy,-^^+h.c)l + ^-^, (73) 

for \ip\ > Mmess, where M^ess is the messenger scale, 7723/2 denotes the gravitino mass, am 
is a constant of order unity and Mp ~ TTig^f^M^^gg with mgoft ~ 1 TeV. 

The lepton number generated through the AD mechanism is estimated as 



— ~ — 75-sm(6t^) 

Tr ^^ flGeVy^ f M \ 



^^^""fiiFtvjtej (10^ j 



1/2 / N 3/2 C^4) 



where Tr is the reheating temperature, 9 is the initial angle of the AD field in the complex 
plane, and I-j/^osI ~ (^^3/2^^)^^^ is the AD field amplitude at the onset of its oscillation. It 
can be checked that thermal effects on the AD field potential is neglected in this parameter 
choice [66, 67]. The AD field fragments into Q-balls after it starts to oscillate, and they 
once dominate the Universe before they decay if < (Tr/3)(|'(/'os|/^p)^5 where is 
the decay temperature of the Q-ball discussed later. In this case, the expression becomes 
ul/s ^ (Trf/4m3/2) sin(66'). Hereafter we regard the lepton asymmetry, denoted by the 
subscript L, as if it is an extra radiation component, which has been denoted by X in 
the previous sections. The neutrino asymmetry in each fiavor is expressed in terms of the 
chemical potential (or the degeneracy parameter) ^y. as 

2^ (75) 



12C(3) \T., 



The neutrino chemical potentials contribute to the extra radiation energy density through 
the relation 

3pL 



i=e,fi,T 



7 \ TT 7 V TT 



(76) 



Note however that the chemical potential of the electron neutrino directly affects the 
helium abundance [68] and hence its contribution to the radiation energy density is con- 
strained as AA^cff ^ 0(0.1) depending on the neutrino mixing angle 6*13 [69]. Thus hereafter 
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we neglect the contribution of the lepton asymmetry to the DR energy density, although 
the effect of its isocurvature perturbation may not be neglected.*^ 

The angular component of the AD field may be light during inflation [63] and it leads 
to the isocurvature perturbation in the lepton asymmetry. The isocurvature perturbation 
of the lepton asymmetry is calculated as 

SL=n cot{n9)S9 - -n'^{S9Y, (77) 
2 

where ^/{SW) = Hia{/{27r\ipi\), with Hiai being the Hubble scale during inflation and n = 6 
in the present model. The DR isocurvature perturbation is then estimated as 



Pdr V 3pDR J Pu \ 3 

' lOSGeV/ VlGeV/ llO^OGeV/ llQi^GeV, 



(78) 



where pdr = Pl + Pu denotes the DR energy density. In the second equality, we have 
neglected the contribution of the lepton asymmetry to the DR energy density for the 
reason discussed above. In the last equality, we have considered only the leading term. In 
the case of Q-ball domination, we obtain 

^--^'^i" (mi^j (^) (loScSv) (lo^itvj • (™) 

Expanding the DR isocurvature perturbation as S'dr — SQb9 + (l/2)5'6i6((5^)^, we obtain 



bf) = — - — racotmt 



See = —n 



1 -^cot\n9) 



Therefore, we obtain the non-linearity parameter for the DR isocurvature perturbation as 

fgpR.SDRgpR _ SgSee S^g t\a'^(u \ rsn 

/nl - + j^,HhL)A,ik,). (81) 

It is easily checked that the first term is of the order of (P^^p/P,^)^ x (ta.n^ {n9) / ANcs) ■ 
Thus the non-linearity parameter can be large enough to be probed for Ps^j^ ~ Pc and 
AiVeff < 1 or tan2(n6') > 1. 

Finally we comment on the Q-ball formation in the present model, which is essential 
for protecting the lepton number from the sphaleron process. After the AD field starts to 
oscillate, the instability develops and Q-balls are formed. We consider the "delayed" -type 



the limit of ANgg <C 1, the total curvature perturbation is conserved for all scales of interest. 
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Q-balls [60], which are formed when the AD field potential becomes dominated by the 
logarithmic term in (73). Then the charge of Q-ball is estimated as [60] 



where /3 = 6 x 10~^, and the radius of the Q-ball is given by Rq ~ '^j^/V Although almost 
all the lepton number created by the AD mechanism is absorbed into Q-balls, they can 
decay into neutrinos from their surfaces. The decay rate is given by Tq ~ Am^^^/ (1927r^(5) 
where A is a surface area of the Q-ball [70]. Then the Q-ball decay temperature, T^, is 
calculated as 

This is well below the electroweak scale. Thus the lepton number liberated by the Q-ball 
decay is not converted into the baryon number.*^ Notice that the Q-ball decay rate and 
hence its decay temperature depends on the charge Q. Thus if the Q depends on the 
initial angle of the AD field 6, the decay rate fiuctuates on large scales and it causes the 
modulated reheating [73] for the case of Q-ball domination. In this case, the magnitude 
of the DR isocurvature perturbation is modified up to an 0{1) numerical factor. In the 
GMSB, however, it is often the case that the ellipticity of the AD field orbit in the complex 
plane is small and the Q does not depend on 9 [60]. Therefore, there is no such an effect. 



6 Summary 

In this paper, we discussed non-Gaussianities in dark radiation isocurvature perturbations. 
Extending our analysis in the previous work [17], we first derived the primordial bispec- 
trum originating from the non-Gaussian isocurvature perturbations in dark radiation. We 
presented primordial bispectra of both the local and quadratic types. As far as primordial 
perturbations have nearly scale-invariant spectra, amplitude of primordial power spectra 
can be parameterized with six non-Gaussian parameters, which consequently measure the 
non-Gaussianities in the mixture of the adiabatic and dark radiation isocurvature modes. 
We also presented CMB bispectrum from these non-Gaussianities, which allows us to fore- 
cast constraints on the non-Gaussian parameters from future CMB surveys including the 
Planck sattelite and a hypothetical CVL survey. While these parameters can be more or 
less constrained from ongoing Planck satellite experiments, there can be still some room 
for future CMB surveys to improve the constraint, and CMBpol [74] and COrE [75] mis- 
sions are desirable to improve the constraints. We referred to SUSY axion models as 

The diffusion process from tlie Q-ball surfaces may transfer the lepton number in the Q-balls into 
surrounding plasma [71, 72] even at the temperature above the electroweak scale. These leptonic charges 
are converted to the baryon number through the sphaleron process, but this amount can be smaller than 
(or comparable to) the observed baryon number [53]. 
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concrete models for non-Gaussian dark radiation isocurvature perturbations and showed 
that they offer distinct signatures on amphtudes in the primordial bispectrum. We have 
also shown that non-vanishing 5'dr, imprinted in the lepton asymmetry, or the neutrino 
density isocurvature perturbation, can be generated through the Affleck-Dine mechanism 
without producing sizable extra radiation energy density. Since observational signatures 
are the same as those in the isocurvature model of the extra radiation component, pri- 
mordial non-Gaussianities in the neutrino density isocurvature perturbation can also be 
constrained by CMB observations. 

Extra radiation with AN^s — 1 will be tested by the ongoing Planck survey with 
high significance and its origin may be identified through the detection of extra radiation 
isocurvature perturbations. Furthermore, isocurvature perturbations in dark radiation 
can offer us unique information for consistent understanding of the early Universe and the 
particle physics theory. 
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